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ABSTRACT 

We present the results of a detailed set of one-zone models th at account fo r the c oupling between 
pulsation and convection following the original prescriptions of iStelliugweril l|l 98fifl . Motivated by 
the arbitrary nature of the input parameters adopted in this theoretical framework, we computed 
several sequences of models that cover a substantial fraction of the parameter space and a longer 
integration time. We found that our models show the same behavior as nonlinear, hydrodynamical 
models, that is, they approach either the limit-cycle stability (pulsational inst ability) or the fixe d 
point (pulsational stability), or present vibrational instability. In agreement with iStellingwerl (|1 flSfljl . 
we find that convection is the main quenching mechanism for pulsational models located across the 
Cepheid instability strip. Moreover, our one-zone models can mimic the pulsational behavior of both 
fundamental and first overtone Cepheids. We also included a turbulent pressure term and found 
that this physical mechanism plays a crucial role in the pulsation characteristics of the models by 
removing the sharp discontinuities along the light and the velocity curves showed by models that do not 
account for turbulent pressure. Finally, we investigated the vibrational and the pulsational stability 
of completely convective models. We consider as the most important finding of the present work the 
identification of a well-defined region in the parameter space where they approach limit-cycle stability. 
The inclusion of turbulent pressure widens this region, thus supporting original suggestions based on 
both linear and nonlinear models of Long Period Variables (LPVs). Several numerical experiments 
performed by adopting different values of the adiabatic exponent and of the shell thickness indicate 
that the coupling between pulsation and convection is the key driving mechanism for LPVs, a finding 
supported by recent theoretical predictions. 

Subject headings: stars: oscillations, variables: Cepheids, Miras 


1. INTRODUCTION 

Variable stars are crucial astrophysical objects since 
they can be used as tracers of st ellar populations in 
the framework of ga laxy evolution (iDolphin et alJl^002l: 
iMonelli et al.l 1200,211 . Moreover, the comparison be¬ 
tween the observables predicted by pulsational mod¬ 
els and the observations themselves s upply an indepen¬ 
dent estimate of the stellar parameters (iBono et al 
iKervella et al.l2001ll . This implies that stellar pulsations 
provide a way to probe regions of the star that would be 
otherwise unaccessible to direct observations. Therefore, 
variable stars provide the unique opportunity to investi¬ 
gate the plausibility of the physical assumptions adopted 
to construct both evolutionary and pu lsational models 
iBono et alJl2002t iKeller fc Woodll2002D . 

Nonadiabatic stellar pulsations are complicated phe¬ 
nomena, even for small amplitudes and purely radial os¬ 
cillations, to which the linear theory is applicable. For 
large amplitudes, the nonlinear effects become impor¬ 
tant and the simplifying assumptions must be treated 
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with caution. During the last decades the work on non¬ 
linear stellar pulsations has been developed along three 
main approaches: simple one-zone models, the formal¬ 
ism of amplitude equations and hydrodynamical mod¬ 
els. Although the third approach provides the most de¬ 
tailed and accurate physical description of the outermost 
stellar layers during the pulsation cycle, it is also true 
that it is difficult to hgure out whether the intrinsic fea¬ 
tures of the models are either a direct consequence of the 
adopted physical assumptions or caused by the numerical 
methods and the spatial resolution adopted to discretize 
the stella r structure into a ser ies of concentric shells — 
see, e.g., iPetroni et alJ ll200,2D and references therein. 



zone models have been introduced with the sole pur¬ 
pose of clarifying the analysis by eliminating possible 
subtle uncertainties introduced by the stability of nu¬ 
merical algorithms and by spatial resolution. This was 
accomplished by considering the stellar envelope as a 
one-zone structure: a single, relatively thin, spherical 
mass-shell concentric with the stellar center. These mod¬ 
els have provided clear understanding of the destabi¬ 
lization mechanisms and of the possible consequences 













































2 


Limit-cycle behavior in convective one-zone models 


of couplings as well as feedbacks between sever al phe¬ 
nomen a associated to ste ll ar va ri ability — se e jBgkeiJ 
ll966|^j_|Usher&^/hitne^_J196^ iSaitou et al.l ill 9 8911 . 
lUnno fc Xiond (llQQ,*!!! . Ild^ et al ] (l1992li . and iTanakal 
< 200111 . 

The simplest one-zone model that accounts for the cou- 
pling between puls ation and convection was suggested by 
iStellingw erJ <198011. Fo l lowing the convective scheme de- 
velooed bv IStellingweril <198211 . he derived a one-zone pul- 
sational model that includes a time-dependent convective 
term. He computed a set of models associated to the 
Cepheid instability strip, integrating them only for a few 
dynamical time scales. The results of this work support 
the generally accepted view that convection is a damp¬ 
ing mechanism, in particular for models located close to 
the cool edge of the instability strip. Although the pulsa- 
tionally unstable region as well as velocity curves were in 
qualitative agreement with empirical data, he also found 
that fully convective models underwent vibrational insta¬ 
bility for certain values of t he parameters. A simi l ar ap- 
proach was als o adop ted bv lUnno fc Xiond <199(1119^ 
and bv iSaitoil <199,811 . but they did not discuss in detail 
the pulsational and dynamical behavior of their models. 

In this work, we present new results and possible ex¬ 
tensio ns of the one-zone model developed bv lStellingwerJ 
<198(111 . investigating the limit-cycle behavior (pulsa¬ 
tional instability) and by accounting for the role of the 
turbulent pressure and of the thickness of the convective 
layer. The paper is organized as follows. In §2 we present 
the one-zone convective model. We place special empha¬ 
sis on the adopted physical and numerical assumptions. 
In §3 we thoroughly invest igate the lim i t-cycle behavior 
of the model suggested bv IStellingweril <198(111 . an issue 
marginally addressed in the original paper. In this sec¬ 
tion, we also focus our attention on the dependence of 
pulsation properties on the shell thickness (§3.1) as well 
as on the turbulent pressure (§3.2) and their role in de¬ 
termining the morphology of light and velocity curves. 
In §4 we discuss the approach to limit-cycle stability of 
completely convective models and briefly describe possi¬ 
ble empirical similarities with LPVs. Finally, in §6 we 
summarize our main results and briefly outline future 
perspectives. In appendix A, we discuss in more detail 
the physical assumptions adopted to include the turbu¬ 
lent pressure in the one-zone model are described. 

2. THE ONE-ZONE CONVECTIVE MODEL 

In the classical theoretical framework of one-zone mod¬ 
els, variable stars have an equilibrium radius Rq and an 
extended shell or envelope of variable radius on top of 
a com pact core of radius Re- As detailed in lStellingwerJ 
ljl98(<l . by accounting for the equation of motion, the en¬ 
ergy equation, and the equation of convective transport, 
the resultant dynamical system describes the evolution of 
the convective upper layer and it has the following form: 

4^ = HX-'^ - x-^ 
dr^ 

4 ^ = - -feX-^u!) ( 1 ) 

dr 

4^ = Cc(W-‘^iLi/2 _ 
dT 

where X and H are the radius and the nonadiabatic pres¬ 


sure normalized to their equilibrium values, while Uc is 
related to the convective velocity and it is defined as 
Uc = U'/U-aAo, where I/mio is the equilibrium mixing- 
length convective velocity. Moreover, the time variable, 
r, is normali zed to the dynami cal (free-fall) time scale. 
As defined in IStellingweril <1 98fifl . the other input param¬ 
eters of the model are: q = mFi —2,d = m{Ti — 2 )/ 2 , 
and c = m — 2 , where Fi is the adiabatic exponent and m 
is the so-called form factor. In the limit of small oscilla¬ 
tions <Stellingwerjll972lll98(ill . the form factor is defined 
as: 


771 = lim 
x^i 


log [{X^ 

log{X) 


3 

1 — 77 ^ ’ 


( 2 ) 


with T] = Rc/Ro- At the same time, b is defined asb = 4+ 
m[n— (s-|-4)(Fi — 1)], where s and n are the temperature 
and the density expon ents in the K ramers’ opacity law 

— for more details see iBakeil <196(111 . 

The main control parameters of the model are C, Cc 
and 7 c, which are defined as follows: i) the nonadiabatic 
parameter — ( — the ratio between the dynamical and 
the thermal time scale; ii) the convective efhciency — 

— the ratio between the dynamical time and the convec¬ 
tive time scale; Hi) the convective/radiative splitting — 
7 c — the ratio between the initial convective luminosity 
and the total initial luminosity, that is, 7 c = Lcq/Lq, 
which implies that 7r = 1 — 7c- 

The reader interested in a detailed discussion con¬ 
cerning the convective time scale in the context of 
one-zo ne models is referred to §4 of iStellingwerfI 
<198611 . More quantitative predictions concerning the 
different timescales connected with turbulent energy 
across t he envelope of var i able stars have be en pro- 
vided by iStelling wer j <1982ll . iBono et <199911 and by 
iFeuchtinger et a,1J <299(111 . 

The values 77 = 1 and s = 3 for the Kramers’ opac¬ 
ity are the standard ones, while the values for Fi, m, 
7 c, C Cc must b e treated with caut ion. The value of 
Fi = 1.1 adopted in IStellingwerfI <198611 is a typical value 
for the 7 -mechanis m operating in the partial ionization 
regions <Coxlll98(lll . The choice of a form factor m equal 
to 10 (77 = ( 3.888) — considered typical of Cepheids in 
iStellingwerj l|1986ll — together with the values chosen for 
the other parameters assure the turbulent stability and 
the sec ular stability for C and Cc in the range [ 0 , 10 ] and 
7c < 1 <Stellingwer Jll986ll . Th e values of these py am- 
eters were arbitrarily chosen bv iStellingwer j <198611 and 
lead to the existence of a pulsationally unstable region 
(“strip”) for values of 7 c < 0.45. From a mathematical 
point of view, we refer to such a behavior as having a 
limit cycle born through a H o pf bif urcation. The sta¬ 
bility analysis in iStellingwerJ <198611 was hampered by 
the fact that individual models were integrated only for 
a small number of dynamical time scales. In particular, 
our calculations show that for t hese cases, n o limit -cycle 
stability exists. Note also that iStellingwer j <198611 ana¬ 
lyzed only the cases for C and Cc < 4, and that an exhaus¬ 
tive study of the dependence of the limit-cycle behavior 
on the input parameters remained to be performed. 

In the current investigation, we have first undertaken a 
thorough study of the dynamics in the (C, Cc) plane for C 
and Cc < 10 and for different values of the convective- 
radiative splitting, 7 c, following the assumption of 
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1.4 



Fig. 1.— Examples of dynamical behavior for 7 c = 0.4. Radius 
time series for stable — (Cc^C) = (3,4) — (^op), limit-cycle — 
= (1.5,1.5) — (middle) and unstable behavior — (Cc^C) = 
(0.5, 0.5) — ('bottom). 

radiative-dominated energy transport ( 7 c < 0.5). We 
have used the same form factor (m = 10) and the same 
initial con dition (Xn.Vg. Hn, Ur„) = (1.4, 0 . 0 , 1 . 0 , 0.7) 
adopted bv lStellingwerJ lll98fill . We have furthermore as¬ 
sumed that a case is pulsationally unstable (limit-cycle 
behavior) when two consecutive maxima in radius vari¬ 
ations differ by less than 10 “^, and pulsationally stable 
(damping oscillations) when the solution asymptotically 
approaches the fixed point (X,V,H,Uc) = (1,0,1,1) 
and its distance to the hxed point becomes smaller than 
10“®. The solution has been considered to be vibra- 
tionally unstable when the rapidly increasing dimension¬ 
less radius amplitude became larger than 15. These 
thresholds are reasonable values. Moreover, we have also 
verified that the resultant qualitative dynamics does not 
depend on the choice of these values. Of course, one 
could use as well the asymptotic perturbation theory to 
compute the properties of the limit cycle. Nevertheless, 
such an approach is beyond the aims of the current in¬ 
vestigation and we consider it as a future line of work. 
For a better understanding of our results, we illustrate 
in Fig. ^ the three types of dynamical behavior we have 
identified. The top panel shows an example of pulsa- 
tional stability, i.e., the initial perturbation decays. The 
bottom panel shows a vibrationally unstable model (the 
initial perturbation grows), while the middle panel shows 
the time series of a case that approaches a limit-cycle 
stability (pulsational instability). Note that for small 
amplitudes (r < 100) the case plotted in the bottom 
panel seems to approach a limit cycle, but for r > 200 it 
becomes clear that the amplitudes are steadily increas¬ 
ing. This finding further strengthens the need for long 
time integrations to assess the type of stability for the 
individual one-zone models. As a final comment on the 
numerical assumptions considered in the present work, 
it is worth mentioning that the dynamical behavior of 
the system is independent of the initial condition, pro¬ 
vided that it belongs to the neighborhood of the fixed 
point (X,V,H,Uc) = (1,0,1,1). Therefore, we did not 
investigate the dependence of the characteristics of the 
time-series on initial conditions. 

3. LIMIT-CYCLE CHARACTERISTICS 

One-zone models can be adopted to perform only 
a qualitative comparison with pulsation properties of 
variable stars. Moreover, this comparison is lim¬ 
ited to observables that one-zone models can account 


for. In particular, the shape of light and velocity 
curves — originally suggested for RR Lyrae stars by 
iStellingwerf fc Donohoel 1 198^ — and the occurr ence of 
hrst overtone pulsators I Stellingwerf et al ] IT^ . The 
most interesting results disclosed by our simulations are 
presented in Fig.| 2 H“e, where the regions of damped os¬ 
cillations are shown in grey, the regions where the ra¬ 
dial displacements approach a limit cycle are marked in 
black, while vibrationally unstable cases are shown in 
white. A glance at the data plotted in Fig. |21 discloses 
that, with the exception of the case in which jc = 0 . 1 , 
an increase in jc causes a decrease in the region of pulsa¬ 
tionally unstable cases and that for 7 c = 0.5 it vanishes. 
These findings sup port the results originally obtained by 
iStellingwerJ 11198^ concerning the damping role of con¬ 
vection. In addition to the radi us (X) and velocity {V) 
curves chosen to be illustrated in lStellingwerJ (1198611 . we 
have added the light (L) curve and in several occasions 
the temperature (T) variations, estimated according to 

L=^ = ^ + 7cA-"t73 (3) 

Lq Lq 

(4) 

40 

where L* denotes the stellar luminosity while the sub¬ 
scripts “r” and “c” refer to the radiative an d convective 
compo nent, respectively, as they result from iStellingwerl 
(119861) . Fig. 0 shows the light curve, (L), the radius 
(A), and velocity {V) time series for selected cases that 
approach limit cycle. Note that to improve the visual¬ 
ization of data plotted in Fig. |31 we subtracted from the 
integration time, r, the time interval spent by individual 
models to approach the limit cycle stability. We have 
chosen to illustrate cases at fixed 7 c and Cc values, and 
C ranging from 0.7 to 6.5, being typical examples of the 
qualitative dynamics encountered. The time series plot¬ 
ted in this hgure show that the cases characterized by 
parameters located close to the edge of the pulsational 
stability region (grey area in Fig.| 2 ) present small ampli¬ 
tudes and sinusoidal changes (panel a). These features 
are quite similar to the behavior that nonlinear, convec¬ 
tive models show close to the blue e dge of the instability 
strip — see, e.g.. lBono et al ] ll2000l) . On the other hand, 
models that are located across the pulsationally unsta¬ 
ble region (panels b, c) show a bump along the rising 
branch of the luminosity curve. Our calculations reveal 
that such a feature is due to the sharp increase in the 
efficiency of the convective motions, close to the phase of 
minimum radius. This finding is also found in nonlinear 
hydrodynamical models. 

One can notice that for cases located close to the vi¬ 
brationally unstable region, the bump becomes very nar¬ 
row (panel d) and pulse-like (panel e). This is a be¬ 
havior that is not observed in actual Cepheids and it 
is not supported by nonlinear, convective models. This 
indicates the limit of the crude physical and numerical 
approximations adopted to construct our models. How¬ 
ever, our analysis suggests that the limit-cycle behavior 
of these one-zone models in the pulsationally unstable re¬ 
gion mimics the behavior of the Classical Ce pheids insta¬ 
bility strip, as argued in the original paper of lStellingweri 
(119861) . Finally, we note that the period of the oscilla¬ 
tions of the cases plotted in Fig. 3 is not unity, even if 
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Limit-cycle behavior in convective one-zone models 



Fig. 2.— The (^c)-plane for several values of 7c: 7c = 0.1 (a)-, 7c = 0.2 (b)-, 7c = 0.3 (c); 7c = 0.4 (d); 7c = 0.5 (e), representing 
pulsationally unstable (black), pulsationally stable (grey) and vibrationally unstable behavior (white) 




Fig. 3. — Limit-cycle characteristics. From top to bottom the panels show: the total luminosity, L; the radiative luminosity, Lr = Li,^/Lo 
and the convective luminosity, Lc = L*c/^o; and the dimensionless radius, X and velocity, V variations. The temporal axis is represented 
as T — ric> where ric is the time when the limit cycle has been reached. The parameters adopted to construct the individual cases were 
chosen from Fig.lJJi: (a) Small amplitude: (CiCc) = (6.5,3); (b) Bump Cepheid: (C>Cc) = (5-5,3); (c) Double-peak Cepheid: (CiCc) = (4,3); 
(d) Steep bump: (C,Cc) = (1.5,3); (e) Pulse-like: (C,Cc) = (0.7,3). 


the dynamical timescale was normalized to unity. This 
str ong nonlinear effec t on the period was already noted 
bv lStellingwerJ lIlQSfiD . and he suggested that it was due 
to a correlation with the pul sation amplitude (see also 
iStellingwerf fc DonohoellT98fiD . 

3.1. The shell thickness 

In this section, we present the results concerning the 
impact of the shell thickness on the existence of a limit- 
cycle behavior and on its characteristics. According to 
classical physical arguments, pulsation al models o nly ac¬ 
count for the envelope of variable stars 100x1198(11) . They 
include a damping, adiabatic region that is typically lo¬ 
cated on top of the nuclear burning region (core), a tran¬ 
sition region, a driving, nonadiabatic region, and the out¬ 
ermost layers (surface). The driving region is the enve¬ 
lope zone where key elements (hydrogen, helium, metals) 


are partially ionized and supply, via the k- and/or the 
7 -mechanism, the pulsation destabilization. In this theo¬ 
retical framework, the base of the envelope in nonlinear, 
hydrodynamical models i s typicallv located a t a few per¬ 
cents of the total radius llPetroni et al.ll20f)^ . 

As far as the one-zone-model approach is concerned, 
different assumptions have been adop ted in the litera¬ 
ture. The one-zone model suggested bv licke et all lll992i) 
for Mira variables accounts for a driving region (piston 
approximation) at A) « 0.2 — 0.4, a transition region 
through which the pressure waves from the interior prop¬ 
agate, and a dissipation region that t hey call mantle. On 
the other hand, the one-zone model of iStellingweril l|1 98(tD 
for variable stars in the Cepheid instability strip accounts 
for a driving convective region located at A) « 0.85 on top 
of a damping, adiabatic region (rigid core). In this the¬ 
oretical framework, the pulsation destabilization is pro- 
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Fig. 4.— The (f, Cc)-plane for several values of 7 c. Eq. Q has been used with rj = 0.75: 7 c = 0.1 (a); 7 c = 0.2 (b); 7 c = 0.3 (c); 7 c = 0.4 
(d); 7 c = 0.5 (e), representing pulsationally unstable (black), pulsationally stable (grey) and vibrationally unstable behavior (white). 



vided by a driving agent (Fi < 4/3). The two differ¬ 
ent assumptions describ e the same e n velop e region if the 
rigid core, assumed bv iStellingwer j l|198fiD is the outer 
boundary of the damping region. Therefore, the shell 
thickness, (1 — r/) in dimensionless formulation, is the 
radial extent of the region located above the boundary 
between the damping (adiabatic) and the driving (nona- 
diabatic) region. 

To investigate the dependence of the dynamic behav¬ 
ior of the system given in Eq. o on the shell thickness, 
we have chosen two values of the shell thickness, one 
smaller and the ot her larger than the value adopted in 
iStellingwerJ lll98(i|l : y = 0.75 and y — 0.92. We illus¬ 
trate in Fig. 0 and Fig. 0 the behaviors in the (C, Cc)- 
plane for jc < 0.5. One can notice that the increase 
in the shell thickness (decrease of 77) leads to the ex¬ 
istence of pulsational instability for cases characterized 
by stronger convection (higher values of 7 c). This ef¬ 
fect is expected, because a decrease in 77 implies an in¬ 
crease in the extent of the driving region. Moreover, 
we have noted that for a fixed value of 7 c and different 
values of 77, the period distribution of the cases that ap¬ 
proach a stable limit cycle peaks at shorter periods as the 
shell thickness decreases. This behavior is also expected, 
since detailed nonlinear, hydrody namical models of RR 
Lvrae llB ono fc S tellingwerall994ll and classical Cepheids 
l|Bono et alJl999D suggest that the regions located below 
the nodal line supply a small contribution to the work 
integral of fir st overtone pulsators. T his supports the 
suggestion bv iStellingwerf et al.l lll987ll to decrease the 
shell thickness of one-zone models to mimic the dynam¬ 
ical behavior of overtone pulsators. 


3.2. The turbulent pressure 

Convection affects pulsation through three factors, 
namely convective energy transfer (thermodynamic cou¬ 
pling), turbulent pressure, and turbulent viscosity (dy¬ 
namic coupling). The effect of turbulent viscosity is to 
convert the kinetic energy of radial motions into ther¬ 
mal energy by means of a turbulent cascade of smaller 
and smaller turbulent eddies. This me ans that the turbu¬ 
lent viscosity is a pure damping facto r l|Stellingweril982l: 
iXiong et alJl998alYecko et alJll998ll . The role that tur¬ 
bulent pressure plays in driving or damping the pulsa¬ 
tion is not straightforward, since the contributions of gas 
and turbulent pressure can not be easi ly separated. Lin- 
ear llYecko et alJll99^ and nonlinear llBono et alJll999li 
convective models of Classical Cepheids suggest that the 
work done by turbulent pressure attains both positive 
(driving) and negative (damping) values in different re¬ 
gions of the envelope. Line ar and nonlinear, co nvective 
LPV models constru cted bv lOstlie & CoxI lll98fiH and by 
iCox fc Ostliel 1 II 993 H indicate that turbulent pressure is 
a driving mechanism, whereas more recent calculations 
bv IXiong et alJ l|1998aji support the evidence that it is a 
damping mechanism. 

The turbulent pressure, Pt, was not include d in the 
original one-zone model bv iStellingwer j l)198(Hl . but it 
was shortly mentioned as being 

A = x-^Ul (5) 

-mo 

where the zero subscript denotes its equilibrium value. 
In this section, we present the results obtained by the 
introduction of the turbulent pressure. The reader inter¬ 
ested in a detailed description of the physical assump- 
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Fig. 7.— Limit-cycle characteristics when the turbulent pressure is introduced. From top to bottom the panels show: the total luminosity, 
L; the radiative luminosity, and the convective luminosity, Lc = L*c/Lo; and the radius, X and the velocity, V variations. 

The temporal axis is represented as r — ric, where ric is the time value when the limit cycle has been achieved. The parameters adopted to 
construct the individual cases are: 77 = 0.888, 7 c = 0.4 and (a.) ((^, (^c) = (4, 0.28); (b) ((^, ^c) = (4, 0.4); (c) ((^, (^c) = (4,1); (d) ((^, ^c) = (4, 4); 
H «,Cc) = (4,9). 


tions adopted to include the turbulent pressure term in 
Eq. O is referred to Appendix A. The new dynamic sys¬ 
tem becomes now: 


^ = (1 - ap)X-9/i + apX-^U^ - A-2 
dr^ 

dH Q:p(r3 — 1) 2d-l dX 2 

— =-m — - A — - 

dr 1 — ckp dr 

- ^ - 1) (6) 

dr 

We have undertaken a parametric study in order to 
investigate the influence of the turbulent pressure on 
the overall dynamics of the system and implicitly on 
the existence of limit cycles. For this purpose, we 


have chosen the same initial condition as in §2, that 
is {Xo,Vo, Ho,Ucg) = (1.4,0.0,1.0,0.7), which led to a 
value of ap « 0.4. The results concerning the plane 
((C, i^'c) are shown in Fig. El using the color code from 
Fig.m The regions characterized by the existence of 
limit cycles for ■jc < 0.5 are more extended than in the 
case without turbulent pressure, as expected because the 
additional pressure adds to the driving mechanism. 

We have also investigated the light and velocity curves 
for the cases characterized by limit-cycle behavior. From 
the hydrodynamical models as well as from plain physi¬ 
cal arguments on the existence of the Cepheid instabil¬ 
ity strip, it is expected that by crossing the instability 
strip from the blue to the red edge, the period increases 
and the amplitude of the oscillations quickly reaches its 
maximum and then slowly decreases. In order to obtain 
a similar evolution of the amplitude from this one-zone 
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model, the limit-cycle region must be bounded in the 
(Cj Cc)-plane by pulsationally stable regions both to the 
left (low values of Cc) and to the right (high values of 
^c)- If bounded by a vibrationally unstable region, the 
oscillations which result from using parameters close to 
this boundary are pulse-like and the amplitude is large, 
as in Fig. Et- One can notice that for the case of zero- 
turbulence pressure, the limit-cycle region is bounded to 
the left by the vibrationally unstable behavior and thus 
the above condition is not satisfied. However, with the 
introduction of the turbulent pressure, there exists a case 
for which the limit-cycle region is bounded in the (C, Cc)- 
plane both to the left and to the right by the pulsation- 
ally stable behavior (Fig. EH)- For this case, we present 
in Fig. [7| the types of light and velocity curves obtained 
in the transition from the hot (blue, low Cc) to the cool 
(red, high Cc) edge, at constant (■ One can notice that 
the amplitude has a peak toward the center of the limit- 
cycle region and then it decreases. Moreover and even 
more importantly, data plotted in Fig. 0 disclose that 
one-zone models that account for turbulent pressure do 
not present, in this region of the parameter space, the 
spurious secondary peak along the rising branch. The 
morphology of light and velocity curves plotted in this 
figure are in qualitative agreement with actual classical 
Cepheids and with nonlinear, convective models. 

4. THE RED VARIABLES 

iStellingwerJ lll98fill remarked a paradox of the present 
one-zone model: although the convection tends to 
stabilize the pulsation in the majority of cases, the 
completely convective case shows instability in every 
criterion. We have verified that indeed the present 
model with fully convective (% = 1) thick shell (m < 8 
or 7f < 0.85) presents vibrational instability for any 
C) Cc < 10. However, for a thin shell, our simulations 
revealed that while no limit-cycle regions exist for the 
cases with 0.5 < 7c < 1-0 and (,(c G [0,10], limit-cycle 
models exist for the case of fully convective shells. 
Fig. EJi shows the (C,Cc)-plane for 7 c = 1.0, ry = 0.888 
and without turbulent pressure. Although the current 
analysis relies on simple one-zone models, this is a very 
interesting finding and we are tempted to attribute 
this region of the plane (C,Cc) to the instability strip 
of LPVs — variable red giants and supergiants — 
which are thought to be significantly nonadiabatic 
and highly convective. There is general agreement 
within the astrophysical community that the intrinsic 
reason for the existence of the red edge of the Cepheid 
instability strip is the damping produced by convection. 
However, theoretical and empirical evidence suggest 
that convection might also be a destabilizi ng mechanism 
in th e cool region of the HR di agram (lFo 2 Lfc_WQQd 
1982t lEdmonds fc Gillilandl Il99(it Ixiong et alJ Il99~ ; 
WoodI 1200(11: IChristensen-Dalsgaard et al.| I200ll: 

Dziembowski et al.l200lUKiss fc Beddinel2f)0^llta et al.1 

200^ When the convective time scale is much longer 

than the dynamical time scale ((c 1), the effect 

of convection is stabilizing, but when the co nvective 
time scale is much shorter, the reverse is true llGou 
11113). Thus, it may provide an important driving 
mechanism in red giants and supergiants. More exactly, 
the convective luminosity in the ionization regions 
of red pulsating variables is expected to exceed 99% 


of the total luminosity, as mentioned by iXiong et alJ 
ill 998ali . They identify a Mira instability strip outside 
the Cepheid instability strip when pulsation-convection 
interaction is taken into account. 

To investigate the effect of the turbulent pressure in the 
case of a completely convective shell, we have identified in 
the (C, Cc)-plane the regions characterized by limit-cycle 
behavior. We illustrate the results in Fig. Et- The region 
of limit cycles for the completely convective case persists 
when the turbulent pressure is introduced. Compared to 
the cases with no turbulent pressure, this region shifts to 
smaller values of the convective efficiency, (c- This shift 
is due to the fact that the minimal perturbation strength 
necessary to drive the pulsational instability is achieved 
at weaker convective driving (smaller (^c), as the rest of 
the driving is now provided by the turbulent pressure. 

Before verifying that convection is the driving agent 
for the completely convective shell, we briefly describe 
the pulsation characteristics of the models constructed 
by adopting jc = 1-0 and compare them with the previ- 


(a) (b) 
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Fig. 8 .— The (^, (^c)-plane for the completely convective shell 
illustrated in the spirit of Fig[^for: jc = 1-0, r) = 0.888 and (a): 
no turbulent pressure, and (b): turbulent pressure included; (c) 
Time series for the luminosity L, radius X, velocity V, pressure 
convective velocity Uc and temperature T for a Cepheid-like 
(solid line: 7 c = 0.2, = 1.2, ^ = 8 ) and a LPV-like variable 

star (dashed line: 7 c = 1 . 0 , Cc = 9? C = 7.5), the latter case was 
extracted from panel (a). 
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Fig. 9.— Variation of the luminosity (L), radius (V), velocity 
(V), pressure (H), convective velocity (Uc), for the case of com¬ 
pletely convective shell close to the vibrationally unstable edge 
= 8.7; = 9). The case of zero-turbulent pressure was consid¬ 

ered. 


ously discussed limit cycles. For comparative purposes, 
Fig.|H|: shows the time series of a “Cepheid” {-jc < 0.5), 
and of a “LPV” model ( 7 c = 1.0). The increase both 
in period and amplitude observed in this figure when 
passing from a “Cepheid”-like to a “LPV”-like supports 
the previous working hypothesis. In this context it is 
worth mentioning that “Cepheid”-like pulsators show the 
well-known phase-lag between the light and the velocity 
curve. The pulsation behavior of the “LPV”-like model 
needs to be discussed in more detail. The data plot¬ 
ted in the top and in the bottom panel show a very 
small phase lag between the light and the temperature 
changes. This evidence ap pears to be suppor t ed by in¬ 
frared spectroscopic data lIHinkle et al ][l98l IT^ of 
pulsating Asymptotic Giant Branch (AGB) stars. How¬ 
ever, the maximum in the light curve takes place before 
the minimum in the velocity curve. This finding is at 
odds with observational data for Mira and Semiregular 
Variables, since the minimum velocity i s not correlated 
with light maxima llH inkle et "mI Il982t iLebzelter et al.l 
l200r)HLebzelter fc Hinklell2002(l . Nevertheless, the com¬ 
parison between theory and observations is hampered by 
secondary periods llHinkle et al.1 
, and indeed in a few objects the 
minimum in the velocity curve takes place later than 
the maximum in t he velocity c urve — y Oph, R Leo 
llHinkle et aLlIlQSll and S Lep llWood et al.1120041 . In 
passing we would also like to draw the attention on the 
correlation between light maxima and convective veloc¬ 
ity, since the maximum takes place just before maximum 
light. We are not aware of macroturbulent velocity mea¬ 
surements in LPVs, and new observations would be very 
useful to constrain the plausibility of the physical as¬ 
sumptions adopted in simple one-zone convective mod¬ 
els. 

For completeness, we illustrate in Fig. the time be¬ 
havior of a model located close to the vibrational in¬ 
stability edge after approaching to limit-cycle stability. 
This example is generic for all the cases with or without 
the turbulent pressure close to the vibrational-instability 
edge for the completely convective shell. The data plot¬ 
ted in this figure show a well-defined bump along the 
decreasing branch of the light curve. This bump be¬ 
comes more pronounced as the values of the parameters 


the occurrence of long 


C and Cc approach the vibrational-instability edge, while 
the peak of the luminosity becomes a pulse-like feature. 
The occurrence of such second ary features is quite typical 
along the light curves of Miras llWood et al.l2003l . More¬ 
over, the occurrence of the luminosity maximum appears 
strongly correlated with the variation of the convective 
velocity {Uc), while the bump is in phase with the change 
of the radial velocity {V). 

Similar to the limit cycles for the weakly convective 
cases associated to the Classical Cepheids, it is highly 
speculative to attribute specific values of the parameters 
C and Cc to different types of LPVs. A more detailed in¬ 
vestigation of the parameter space is mandatory to pin¬ 
point the pulsation properties of completely convective 
one-zone models that mimic an “RGB”-like or a “Mira”- 
like behavior. 


4.1. The 'y-mechanism 

As a natural further step into clarifying the origin of 
the pulsationally unstable cases for the completely con¬ 
vective shell, we have investigated the role of the adia¬ 
batic exponent, Fi, as driving agent. From the math¬ 
ematical perspective of the model, the question resides 
into establishing the relative effect on the system dynam¬ 
ics of the parameter 7 c, on one hand, and of the parame¬ 
ter Fi, on the other hand, included in the definition of the 
parameters q = mTi — 2 and d = to(Fi — 2)/2. The case 
of completely convective shell translates for this one-zone 
model into the disappearance of the At-mechanism, with 
the driving being supplied only by the 7 -mechanism. To 
test whether the pulsational instability of our convective 
models is due to the 7 -mechanism or to a convection- 
induced driving mechanism, we studied the evolution of 
the limit-cycle region in the (C, Cc)-plane as Fi is varied. 
The results obtained so far have used the value Fi = 1.1 



Fig. 10.— The (i^, Cc)-plane for the completely convective shell 
( 7 c = 1.0) and different values of rj and Fi giving birth to several 
behaviors: limit cycle {filled circles), pulsational stability (dots) 
and vibrational instability {white regions), (a) rj = 0.92, Fi = 
1.1; (b) I? = 0.888, Fi = 1.1; (c) q = 0.888, Fi = 1.2; (d) q = 0.888, 
Fi = 1.3. Eq. □ has been used. 
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as the one employed bv iStellingwerJ lll98fiD . We conjec¬ 
ture that if the obtained cases of self-sustained oscilla¬ 
tions are a product of convection-induced driving, then 
the 7 -mechanism has very little, if any, influence on their 
properties. Fig. nnii and b display two cases of different 
shell thickness concluding, as in §3.1, that a thicker shell 
implies a wider range of parameters leading to pulsa- 
tional instability. In the completely convective case, this 
increase causes a shift of the limit-cycle region to lower 
^c- The same shift occurs by increasing the Fi value at 
constant shell thickness, or, in other words, by decreas¬ 
ing the efficiency of the 7 -mechanism iFig.ll Oh.c.rl). This 
means that convection and turbulent pressure are driv¬ 
ing mechanisms for fully convective models, since the re¬ 
gion where these models display a limit cycle stability 
marginally depends on the adopted Fi valu e. This hnd- 
ing sup ports the theoretical predictions bv iXiong et al.1 
Il1998ati concerning the driving provided by the coupling 
between pulsation and convection, but we also find that 
turbulent pressure is a destabilizing mechanism instead 
of a damping factor for pulsation. However, current mod¬ 
els assume that the physical properties of the driving re¬ 
gion can be described by a Kramers’ opacity law and 
by a very crude equation of state. Therefore, there is 
no guarantee that current models properly account for 
the thermodynamical and dynamical properties of the 
different partial ionization zones in low density envelope 
regions. 

5. CONCLUSIONS 

In this paper we have thoroug hly analyzed t he on e-zone 
convective model introduced bv iStellingweril lll98(ill . The 
model appears in the form of a system of 4 ordinary dif¬ 
ferential equations where the variables are the radius of 
the shell, the velocity, the pressure and the convective ve¬ 
locity. The nonadiabaticity resides in the pressure which 
is considered as a nonadiabatic perturbation of the refer¬ 
ence pressure and, thus, can be considered as a nonadia¬ 
batic variable of the system. The model accounts for the 
self-excited oscillations by adopting a value for the adia¬ 
batic exponent, Fi, close to unity. The main parameters 
of the system are the fraction 7 c of convective luminosity 
with respect to the total luminosity, and the time-scales 
ratios, and Cc, which are a measure of the dynamical 
time scale to the thermal time scale and to the convective 
time scale, respectively. We have extended the model by 
considering explicitly the role of the turbulent pressure 
and the case of completely convective shell. 

The one-zone models constructed by iStellingwerJ 
<ll98(itl were integrated only for a few dynamical time 
scales and, therefore, the approach to limit-cycle stabil¬ 
ity remained still to be investigated. To properly iden¬ 
tify the region of the parameter space which shows a 
limit-cycle behavior (pulsational instability), we have se¬ 
lected a parametric space given by 7 c G [0,1), C ^-nd 
Cc € [0,10]. Our parametric study revealed well-defined 
regions where limit cycles exist, born through the Hopf 
bifurcation. For a typical shell thickness of 11%, pulsa¬ 
tional instability was encountered only for the radiation- 
dominated cases ( 7 c < 0.5), while the increase in the 
shell thickness displaces this limit to upper values of 7 c. 
Th e existence of thes e upper values supports the work 
bv iStellingwer j lll98fill and further strengthens the role 
played by convection as a damping mechanism. More¬ 


over, we found that the shell thickness is, as expected, 
positively correlated with the period distribution of pul- 
sationally unstable cases. Interestingly enough, models 
constructed by adopting thick and thin shells obey differ¬ 
ent PL relations that mimic the behavior of fundamental 
and first overtone Cepheids. 

Additionally, we have undertaken a parametric study 
to investigate the influence of the turbulent pressure on 
the overall dynamics of the system and implicitly on the 
existence of limit cycles. The turbulent pressure appears 
to be a driving mechanism as the regions where the mod¬ 
els approach a limit-cycle stability are more extended 
than for models that do not account for the turbulent 
pressure. Moreover and even more importantly, we find 
that one-zone models that account for turbulent pressure 
do not show a spurious peak along the rising branch of 
light curves. The inclusion of this term is also supported 
by the occurrence in the (C, (Cc)-plane of a pulsationally 
unstable region bounded to the left and to the right by 
two pulsationally stable regions, as expected from phys¬ 
ical considerations and hydrodynamical models. More¬ 
over, the models located in this region display a mor¬ 
phology of the light curves and the velocity time series 
quite similar to those of nonlinear, hydrodynamical mod¬ 
els and actual Cepheids. 

As a natural continuation of the work bv lStehingwerfl 
(|1 98(ili . we have investigated the vibrational and pulsa¬ 
tional stability of completely convective models. Our 
predictions disclose a well-defined region of the param¬ 
eter space where these red models approach a limit- 
cycle stability. The turbulent pressure appears to be a 
driving mechanism, a hnding that supports the r esults 
original brought fo rward by lOstlie fc CoxI l|198(ill and 
iCox fc OstlielTrflQ.lIl on the basis of both linear and non¬ 
linear LPV models, but at odds with predi ctions based 
on line ar, convective models provided by IXiong et alJ 
l|1998a|i . We computed several sequences of models by 
adopting different values of the adiabatic exponent and 
of the shell thickn ess. We found in agreement with 
IXiong et a, 1 .l (|1 998a|i that the coupling between pulsation 
and convection is the key physical mechanism that drives 
the pulsation instability in these simple structures. We 
have also performed a qualitative comparison with empir¬ 
ical properties of LPVs. We have found some similarities, 
but only one feature partially agrees with empirical data: 
the maximum in the light curve is not correlated with the 
minimum in the velocity curve. However, the comparison 
with observations might be hampered by the occurrence 
of long secondary periods. We are not aware of detailed 
measurements of macroturbulent velocity along the pul¬ 
sation cycle of LPVs, however our models suggest that 
the maximum in the light curve takes place soon after 
the maximum in the convective velocity. 

Current one-zone models account for the coupling be¬ 
tween convection and pulsation and for turbulent pres¬ 
sure. We investigated the sensitivity to free parameters 
and to the physical assumptions adopted to construct 
the models. However, the treatment and the inclusion of 
these physical ingredients rely on crude physical approxi¬ 
mations. A more detailed investigation is required before 
we can assess whether our current theoretical framework 
might mimic the complex pulsation behavior of LPVs 
and this implies the use of the asymptotic perturbation 
theory to compute analytically the properties of the limit 
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cycle. 

The pulsational behavior disclosed by these simple 
models for fully convective models will be addressed in 
a forthcoming paper on the basis of nonlinear, hydrody- 
namical models (Munteanu et al. 2005, in preparation). 

It is a pleasure to thank P. Wood for a detailed reading 
of a draft of this manuscript and for several suggestions 
concerning the comparison between theory and observa¬ 
tions. This work has been partially supported by the 

APPENDIX 

THE TURBULENT PRESSURE 

In order to include the turbulent pressure and its associated energy terms in the present one-z one model, one mus t 
start from the equation of momentum and energy conservation for the treatment of convection of iStellingweri lll982ll : 

(Al) 

_iv(F, + Fe + Ft). (A2) 

P 

Here a s well as throughout our work, we have used the notations employed in the original investigations. In iStellingweri 
11982^1. a ny quantity was written as a; = {x) + x', where (x) was the mean quantity and x', the fluctuating part. In 
En. IIA2ll . D/Dt = {d/dt + (u) ■ V) is the Lagrangian time derivative with u being the convective velocity, E is the 
specific internal energy, U = l/p is the specific volume, P is the thermodynamic pressure and 


= __ v(P + Pt)-V<I> 

Ut p 


A(e + e.) + (p + 


= 

Dt 


MCYT grant AYA2002-04094-C03-01, by the European 
Union FEDER funds, and by the CIRIT. One of us, G.B., 
acknowledges for partial support INAF within the frame¬ 
work of the project: “The Large Magellanic Cloud a lab¬ 
oratory for stellar astrophysics “, and PRIN 2003 within 
the framework of the project: “Continuity and Discon¬ 
tinuity in the Galaxy Formation“. A.M. acknowledges 
a scholarship of UPC which allowed her to conclude the 
present work. 


Et^\{{u'f). 

(A3) 

Pt^p{{u'f), 

(A4) 

Ee = pGp((w'T0), 

(A5) 

Ft ^ 

(A6) 


represent, respectiv ely, the conv e ctive energy, the turbulent pressure, the conyectiye and turbulent kinetic fluxes, as 
they are dehned in iStellingwerfl l|1 982fl . These equations must be adapted to the one-zone model. The momentum 
equation translates into 


(fX 


= (I - QfpjA-^/i -b apA-^U<f - A 


--2 


^2 — IL T 

where Eq. © has been used and the parameters are q = mTi — 2, c = m — 2 and 


Pi 


to 


^0"’"C^c^o 


Pn -b Ptr 


0 -I- vi-o 

To ease the calculation, the energy equation can be diyided into 




(A7) 


(AS) 


d{Lr + Lc) _ TiP dp I dP 

dm — 1) dt p(r 3 — I) dt 

dm p^ dt dt 

In the determination of the conyectiye and turbulent luminosities, one can use the conseryatiye choice of of Eq. (27) 
of IStellingwerfl l)I98(Hl and Eq. IjAbp . The calculations lead to an expression for Lt identical to that of Lc, that is 
Lt = X-’^U^, while Et = Uf/2. Using Eqs. <A9HAinil . one can get the final equation of energy conseryation taking 
into account the turbulence: 


(A9) 

(AIO) 


dh 

dr 


p(r3 -1) 

I Olr\ 


X 


2d-l 


dX 

dr 




-X'^^U'^) - {l,X% 


rblS + i 


- Cc {X<^h^/^u, - 

G) 

+ (I-7r)A-^C/3-I), 


(All) 
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where d = m(Ti — l)/2 and 7 ]. = All the other symbols used in Eq. irmi have their usual meaning. The 

equation for the convective velocity coincides with the zero-turbulence case as the approximation made for its recovery 
uses of the temperature as a function of the thermodynamic pressure only. In the case in which the turbulent energy, 
Et, is neglected in Eq. EH), Eq. lUlTll becomes 


dh 

dr 


= — m 


ap(r3 -1) 

1 Oo 


A 


2d-l 


dX 

dr 


Ui - 


- CX^‘‘ -h (1 - 7r)A-=[73 _ 


(A12) 


Thus, to a first approximation, the one-zone convective model with turbulent pressure is described by Eqs. ITtIi and 
( 1112 ), while for the convective velocity Eq. o remains valid. 
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